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First principles calculations of the E5(310)[001] symmetric tilt grain boundary in Cu with Bi, Na, 
and Ag substitutional impurities provide evidence that in the phenomenon of Bi embrittlement of 
Cu grain boundaries electronic effects do not play a major role; on the contrary, the embrittlement 
is mostly a structural or "size" effect. Na is predicted to be nearly as good an embrittler as Bi, 
whereas Ag does not embrittle the boundary in agreement with experiment. While we reject the 
prevailing view that "electronic" effects {i.e., charge transfer) are responsible for embrittlement, 
we do not exclude the role of chemistry. However numerical results show a striking equivalence 
between the alkali metal Na and the semi metal Bi, small differences being accounted for by their 
contrasting "size" and "softness" (defined here) . In order to separate structural and chemical effects 
unambiguously if not uniquely, we model the embrittlement process by taking the system of grain 
boundary and free surfaces through a sequence of precisely defined gedanken processes; each of 
these representing a putative mechanism. We thereby identify three mechanisms of embrittlement by 
substitutional impurities, two of which survive in the case of embrittlement or cohesion enhancement 
by interstitials. Two of the three are purely structural and the third contains both structural and 
chemical elements that by their very nature cannot be further unravelled. We are able to take 
the systems we study through each of these stages by explicit computer simulations and assess the 
contribution of each to the nett reduction in intergranular cohesion. The conclusion we reach is that 
embrittlement by both Bi and Na is almost exclusively structural in origin; that is, the embrittlement 
is a size effect. 

PACS numbers: 61.72.Bb, 05.70.Np, 61. 72. Mm, 68.35.Dv 



Introduction 

Impurity induced embrittlement of metals is an age 
old phenomenon. For the case of bismuth in copper, it 
was apparently first recorded in the scientific literature 
in the 1870sfi becoming a matter of both manufacturing 
importance and scientific curiosity in the 20th century. 
As far as we know the suggestion that segregation to 
grain boundaries was the culprit was first made by Cyril 
S. Smith in a comment to a paper in 1929^ and this expla- 
nation is now generally accepted as the truth. Perhaps 
the most dramatic event resulting from grain boundary 
segregation embrittlement was the failure of the Hink- 
ley Point 'A' nuclear reactor in 1969.^'^ Given that fast, 
brittle failure is a result of the segregation of insoluble 
impurities to grain boundaries, the physical, microscopic 
origin of the loss of cohesion is still a matter of intense 
debate and controversy. And origins may indeed be dif- 
ferent for different metal and impurity combinations. At 
issue is the question, how can impurities segregated to 
grain boundaries at concentrations of a monolayer (ML) 
or less result in reduction in fracture toughness amount- 
ing to a ductile to brittle transition? 

Opinions are currently divided into two schools of 
thought depending on whether it is believed that changes 
in the chemical bonding between host atoms due to the 
presence of impurities in their neighbourhood is the cause 
of the embrittlement. The alternative to such electronic 
mechanisms can loosely be described as a size effectmi^ 
The electronic mechanisms of grain boundary embrittle- 



ment can in turn be divided into grain boundary weak- 
ening and grain boundary stiffening arguments. In the 
former, the impurities withdraw some electronic charge 
from neighbouring host atoms which makes their cohe- 
sion across the boundary less strong.* In the latter, the 
impurity presence at the boundary leads to the formation 
of strong directional bonds which make it more difficult 
for dislocations to glide during decohesion. This leads to 
local stress piling up in the vicinity of the grain boundary, 
and the crystal falls apart in a brittle manner i2iiSiii 

The other group of mechanisms, in which the effects 
of electronic charge redistribution do not play a de- 
cisive role, has been relatively little studied. 12, As we 
have recently argued, the well known effect of the em- 
brittlement of grain boundaries in copper by bismuth 
is a perfect illustration of such embrittlement. Very 
briefly our arguments were supported by three indepen- 
dent pieces of evidence obtained during a study of the 
S19a(331)[110] symmetric tilt grain boundary with seg- 
regated Bi. Firstly, the electronic structure showed no 
evidence of the "white line" in the electron energy loss 
near edge structure (ELNES) reported experimentally to 
support the charge transfer principle. (The existence of 
the white line has also been called into question on ex- 
perimental grounds.^'') Secondly, elastic constants of bulk 
Cu with impurities (Bi or Pb) strongly depend on vol- 
ume, but not on the presence of the impurity. Thirdly, a 
Kanzaki analysis is consistent with a central force model 
of Cu-Bi bonding, in which case the ease of dislocation 
emission from a crack tip should not be affected by the 
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presence of Bi. 

An opposite view was simultaneously expressed by 
Duscher et al.^ according to whom the role of Bi at a 
grain boundary is to donate some of its electronic charge 
to the neighbouring copper atoms rendering them "zinc 
like" and resulting in the weakening of the Cu-Cu bonds 
across the boundary. This point of view then falls into 
the class of electronic mechanisms. Their conclusion was 
based on a study of the E5(310)[001] symmetric tilt grain 
boundary with an areal density of Bi atoms of 2.9 nm^^ 
(approximately 0.25 ML of Bi at the grain boundary 
plane). Evidence of charge transfer between Bi and Cu 
taking place was based upon experimental ELNES spec- 
tra. In addition, the authors of ref. mention that the 
calculated Bi induced expansion of the grain boundary 
is far too small to be the reason for the grain boundary 
becoming weak if a "size effect" argument were correct. 
On the other hand, our recent arguments to support a 
size effect included the finding of a large grain boundary 
expansion at the E19a Bi segregated grain boundarypi^ 

We decided to address the controversy further by turn- 
ing to the the S5 copper grain boundary so as to make 
a direct comparison with the conclusions of Duscher et 
ali^ Our theoretical calculations presented for the first 
time here, arrive at a very different conclusion. The 
grain boundary expansion turns out to be larger in our 
study and we contend that the grain boundary studied 
in ref. was not embrittled — indeed these authors did 
not test this point either experimentally or theoretically. 
We find that at least one monolayer of Bi at this grain 
boundary is needed for embrittlement. Under those cir- 
cumstances the S5 grain boundary is particularly suscep- 
tible to Bi embrittlement.^^ Comparison of the effect of 
Bi with that of Na and Ag provide further support for the 
principle that in the weakening of Cu grain boundaries 
by Bi, charge redistribution plays a negligible role. 

More importantly, in this paper we suggest a unified 
approach within which the effects of an impurity can be 
split into those related to mechanical distortion and to 
chemical identities of impurity and host atoms. We then 
proceed with our ab initio data to demonstrate how these 
combine together to explain the deterioration of grain 
boundary cohesion for purely structural reasons. 

The structure of this paper is as follows. We briefiy 
summarise the connection between fracture mechanics 
and thermodynamics of decohesion in section in or- 
der to make clear how the quantities we calculate re- 
late to fracture toughness. Section^offers a general ap- 
proach as to how the various contributions which weaken 
or strengthen a given boundary can be quantified using 
explicit calculation of the changes in total energy during 
the stages of a thought experiment. Technical details of 
the calculations, parameters used and tests of the method 
are described in section IIIII We present our results for 
the pure Cu grain boundary and grain boundaries seg- 
regated with Bi, Na and Ag in sections IIVI Ivl and IVII 
Section fVIII aims to identify the dominant mechanism re- 
sponsible for grain boundary embrittlement by Bi and 



Na in terms of size, stiffness and chemical properties of 
the impurity atoms. The section follows the method out- 
lined in section ^ Our conclusions may be found in 
section rvml 

I. THERMODYNAMICS OF GRAIN 
BOUNDARY EMBRITTLEMENT AND 
RELATION TO TOTAL ENERGY IN 
ELECTRONIC STRUCTURE CALCULATIONS 

Grain boundary decohesion is a complex process in- 
volving propagation of a crack under load and possibly 
simultaneous emission of lattice or partial dislocations. 
Its proper description requires atomistic simulations at 
the mesoscopic level which is out of the reach of ab ini- 
tio techniques. The approach that makes the problem 
manageable is due to Rice, Thomson and Wang>i2ii2i An 
atomically sharp crack at a grain boundary will advance 
either in brittle or ductile fashion depending on whether 
the minimal energy per unit area of crack advance (en- 
ergy release rate) required to emit a single dislocation, 
^disb is higher or lower than the energy release rate as- 
sociated with brittle cleavage of the crystal along the 
grain boundary, C/cioav'^^^ If Scicav < ^disi, the crack re- 
mains atomically sharp during decohesion which is there- 
fore brittle, otherwise the system begins to emit disloca- 
tions which blunt the crack. In practice, ^disi can be ac- 
cessed through the anisotropic theory of crystal defects; 
for typical, pure Cu grai n boundaries, ^disi falls in the 
region 1-2 J/m^ (ref. |23|). Our task is thereby reduced 
to determining whether segregation may reduce Gcio&v to 
a value significantly lower that this, hence triggering a 
ductile to brittle transition. 

As regards Scioav , in the limit of perfectly elastic brittle 
fracture, neglecting lattice trapping and energy dissipa- 
tion through phonon emission, it represents the reversible 
work needed to convert unit area of grain boundary into 
twice that area of free surface. The analysis is greatly 
complicated by impurity segregation compared to frac- 
ture in pure materials. This is because, whereas it may 
be taken that the impurity at the grain boundary prior to 
fracture is in equilibrium with that in the adjoining crys- 
tals, having an excess concentration Fgt per unit area, 
the concentration profile at the surface cannot generally 
be predicted a priori. As pointed out by Seah,^ and by 
Hirth and Rice^l one may imagine two limiting cases. In 
the first we suppose that the growth of the crack is slow 
compared to the rate of concentration equilibration be- 
tween the surfaces produced and their underlying bulk. 
In that case the grain boundary and surface excess con- 
centrations of impurity are at all times in accord with 
the equilibrium adsorption isotherms and depend only on 
the temperature and bulk concentrations, and wc identify 
Gcicav with the work of adhesion^ 

Wm = 27, - Jgb. (1) 

Here, jgi, and 7, are the energies per unit area of the 
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grain boundary and surface whose impurity concentra- 
tions are such as to be in equihbrium with the bulk. 
Equation is known as the Young-Dupre equationSi 
when used to analyse the equilibrium contact angle at a 
droplet on a surface. It is relevant in fracture only if the 
impurity diffuses very rapidly in the bulk, for example 
in instances of hydrogen embrittlement. The second case 
is applicable here, namely that the growth of the crack 
is so rapid that the impurity remaining at the surfaces 
has no time to come into equilibrium with the underly- 
ing bulk. Indeed at best some surface diffusion may oc- 
cur as the crack is opening and one may suppose a local 
equilibriurr& is achieved in the surface layer of host and 
impurity atoms^ In that case one identifies Qdcav with 
the reversible work of separatior?- under the constraint 
that bulk diffusion is suppressed. As now the number 
of particles, and not their chemical potentials, are fixed 
during the separation we obtain 

where A is the surface area, G is the excess Gibbs free 
energy of a representative piece of material containing 
either the grain boundary (before decohesion) or the sur- 
face (after decohesion), r^f, is the impurity excess at the 
grain boundary before cleavage, and we assume that the 
impurity is evenly distributed between the two newly cre- 
ated surfaces. The excess in Eq. (|2J is taken with respect 
to the adjoining bulk phases. 

To express equation ^ in terms of surface and grain 
boundary energies one can follow the analysis of Hirth 
and RiceiSi Under the condition that the bulk impurity 
concentration is dilute we may write 

w^Mb) = 27,(irg,) - + - tii)rgb. (3) 

This may be useful in revealing two differences between 
the "slow" and "fast" limiting cases. In finding equa- 
tion Q it is necessary to regard the surface atomic layer 
as a thermodynamic system, uncoupled from the under- 
lying bulk. The impurity atoms have a chemical poten- 
tial /i,^, different from their chemical potential /Xj in the 
bulk and at the grain boundary prior to fracture. One 
can then interpret ^s{\rgi,) as the surface energy of a 
crystal whose bulk composition is such that its equilib- 
rium surface excess is ^Tgi,. Secondly, the final term in 
equation © may be interpreted as the energy released to 
the impurity atoms as a result of the change in chemical 
potential on finding themselves out of equilibrium with 
the underlying bulk. This illustrates the point that one 
would always expect the work of adhesion to be smaller 
than the work of separation. 

Both of the limiting cases embodied in equations 
and (PI , may alternatively be formulated in terms of the 
work of adhesion or separation of the equivalent grain 
boundary in pure material A, say, Wsop(A), plus terms 
accounting for the presence of impurity.^'®. In the case 
of "fast" fracture, if one neglects a small contribution to 



Wsop associated with the change of configurational en- 
tropy when a grain boundary with impurity splits into 
surfaces, equation (O becomes 

eclcav = Wsep(B) = W^,ep(A) -G,eg(Bs) +Gscg(B) , (4) 

where B denotes material with segregant, as in figure J^l, 
and Gsog(Bs) and Gseg(B) are the segregation free ener- 
gies per unit area of surface and grain boundary respec- 
tively. Equation (0J) may be interpreted to imply that an 
impurity which segregates more readily to a free surface 
than to a grain boundary will cause a reduction in work 
of separation, or conversely will lead to grain boundary 
cohesion enhancement^ It is worth repeating that a re- 
duction in work of separation does not necessarily imply 
embrittlement which in the simplified picture advocated 
here requires that ^cioav < 0disi- 

In our ah initio supercell calculations the effect of tem- 
perature is neglected hence we approximate the Gibbs 
free energies Ggb and Gg in jSJ with total energies Egb 
and Eg of the supercells containing the grain boundary 
and the surface, respectively. Correspondingly, the seg- 
regation free energies in Q are estimated as the energy 
required to remove all the impurity from the interface 
and distribute it in the bulk. 

The assumption that impurity atoms distribute evenly 
between the surfaces in equation (0) may be removed in 
a straightforward way. Modification to any arbitrary dis- 
tribution is allowed as long as the surface excesses sum up 
to that at the grain boundary. For example, if one of the 
surfaces stays bare, 2Gs{^rgb) in equation (0) must be 
replaced with GsiEgb) -t- Gs(0). In principle, one should 
choose the distribution that minimises the sum of the 
"surface" Gibbs free energies Gg and therefore the work 
of separation (0) . In our study, we consider two limiting 
cases having either equal distribution, or all impurities 
residing on one of the two surfaces. The difference is gen- 
erally rather modest; experiment usually shows an even 
distribution of Bi between the surfacesi^ 



II. STRUCTURAL AND CHEMICAL 
EMBRITTLEMENT 

In this section we define concepts of "structural embrit- 
tlement" and "chemical embrittlement" precisely. While 
there is no unique definition that can unambiguously sep- 
arate structural from chemical effects, we seek a defini- 
tion that is practical and could be used to rationalise 
results for many different systems. For this reason a gen- 
eral definition should be unambiguous and independent 
of the components, including both interstitial and sub- 
stitutional impurities. For future applications we formu- 
late such a general definition, although the calculations 
described in this paper are for substitutional impurities 
only. Furthermore, some modification to the following 
discussion will need to be made to describe the "slow" 
fracture of equation 
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FIG. 1: General scheme for discussion of embrittlement mech- 
anisms due to boundary weakening by impurity segregation. 
A and B denote the fully relaxed equilibrium unsegregated 
and segregated grain boundaries, respectively. Boundary C is 
created by substituting impurity atoms in B (black spheres) 
with vacancies (white spheres) without further relaxation, and 
D is created from A in the same way. As— Ds denote the 
free surfaces into which grain boundaries A— D cleave. Path 
A-^D^C^B refers to substitutional impurities, whereas 
path A— >C^B applies to interstitial impurities. 



A. Intermediate configurations 

In order to proceed we need to define a number of 
hypothetical reference systems (see figure QJ, which we 
denote A, B, C and D. Each system is in the form of a 
grain boundary of unit area and separated surfaces, the 
latter being distinguished from the grain boundary by 
the suffix s. 

These systems refer to a given host and impurity 
combination, a given grain boundary and a given dis- 
tribution of impurity atoms at the boundary. A system 
comprising separated surfaces should have exactly the 
same total numbers of atoms of each species as the 
corresponding grain boundary, so that bulk quantities 
will cancel out of all the relevant energy differences, as 
described by works of separation. To avoid a complicated 



treatment of ensembles that would be required to deal 
with more concentrated solutions, we assume that on 
the length scale of a supercellS^ the bulk can be treated 
as pure host. The systems are defined as follows: 

A: The relaxed grain boundary of the pure host. 
Ag! The relaxed, separated surfaces corresponding to A. 
B: The relaxed grain boundary of the host with segre- 
gated impurity atoms. 

Bg! The relaxed, separated surfaces corresponding to B. 
C: The same atomic structure as B, with impurity 
atoms removed from the impurity sites defined in B but 
no further relaxation. We refer to the imprint of these 
absent impurities as "ghosts" since they define forces 
that maintain the host atoms in a distorted structure, 
but make no other direct contribution to the energy. 
For substitutional impurities a ghost is a more or less 
strained vacancy. For interstitial impurities a ghost is 
a centre of strain in the host lattice without missing 
atoms. 

Cs- The same atomic structure as Bg, with impurity 
atoms removed from the impurity sites defined in Bg 
but no further relaxation. 

D: This system is only relevant for substitutional 
impurities such as those discussed in the rest of this 
paper. It is the same structure as A, with host atoms 
removed from the impurity sites defined in B but no 
further relaxation. 

Dg! The same atomic structure as Ag, with host atoms 
removed from the impurity sites defined in Bg but no 
further relaxation. 

Some points of clarification may help to visualise these 
systems. Note first that A, Ag, B and Bg label fully re- 
laxed structures, whereas C, Cg, D and Dg label notional 
structures we set up in the computer using the atomic 
positions found for the relaxed systems. The works of 
separation for each system are defined as the differences 
in total energy, for example for system A: 

WscpiA) = E{As) - E{A) . (5) 

Although wc would prefer E to mean the Gibbs free en- 
ergy, for simplicity of calculation we use as a surrogate 
the total energy at OK (see section IJ). A more exact cal- 
culation could be made in principle, to include phonon 
free energies in the harmonic approximation. 

If the impurity is detrimental to cohesion, the work of 
separation of our embrittled boundary, B, is less than 
the work of separation of the clean boundary, W4ep(A) — 
M4op(B) > 0. As summarised in section^ this picture of 
grain boundary weakening is equivalent to the statement 
that surface segregation is energetically preferred to grain 
boundary segregation, that is, from equation I^J, 

Wsop(A) - WsapiB) = Gseg(B,,) - G,eg(B) > . (6) 

These formulations of the difference show how the bulk 
energies, including the energy contribution of dissolved 
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impurities, cancel out. We are now in a position to dis- 
cuss what we mean by a "size effect" or "structural effect" 
as opposed to a "chemical effect" . 

B. Size effects and chemical effects 

Leaving aside for the moment chemical effects, we 
might regard a misfitting impurity in the bulk as stor- 
ing energy mechanically, like a compressed or stretched 
spring. The stored energy is fully released when the im- 
purity segregates to a surface, but only partially released 
when the impurity segregates to a grain boundary. This 
would be a simple mechanical picture of the origin of 
the embrittlement. However, it is unsatisfactory without 
specifying what happens to the stored energy. Consider 
the case of an oversized and therefore substitutional im- 
purity. If it found a sufficiently loose site at the grain 
boundary, it could eliminate the strain energy and lower 
the grain boundary energy by segregating there. In this 
case there would be no further reduction in stored energy 
on cleavage. The net effect would therefore be an in- 
crease in Wscp- While the stored energy in an oversized, 
compressed impurity is likely to be a driving force for 
getting it to a boundary, what matters for embrittlement 
is how much of this energy is stored in the boundary, to 
be released on cleavage. It is stored in three ways: (i) 
as broken host-host bonds, because the impurity sub- 
stitutes for the host; (ii) as stretched host-host bonds; 
and (iii) in the compressed or expanded impurity atoms 
themselves. 

If the impurity is interstitial, the first of these is irrel- 
evant, since we suppose that no host-host bonds in the 
boundary are broken by the presence of the impurity. 
The first two contributions are illustrated in figure ^ by 
the paths labelled "Host Removal (HR) Substitutional 
Structure (SS)" and "Interstitial Structure (IS)" respec- 
tively. System D was introduced solely for the purpose 
of quantifying the magnitude of the substitutional effect, 
by removing the effects of bond-stretching. Accordingly 
we define the following contributions to the path from A 
to C for substitutional structures: 

HR = Ty,ep(D)-iy,ep(A) (7) 

SS = T^sep(C)-W;ep(D), (8) 

and for interstitial impurities 

The final step around the diagram in figure ^ labelled 
"Chemical + Compressed Impurity (CC)" takes the sys- 
tem from C to B: 

CC = Wscp(B)-W^,ep(C). (9) 

This is where the chemistry of the impurity comes into 
play. Notice that in this step the host atoms are in the 



same positions with and without the impurity, so no fur- 
ther energy change is associated with the energy stored 
in the host lattice. However, if we are thinking mechan- 
ically, there may still be some energy released by vol- 
ume relaxation of the impurities in going from the grain 
boundaries to the surfaces. We might expect this to be 
more significant when the grain boundary concentration 
is low, in which case there are still some host-host bonds 
holding the impurities in a state of compression, com- 
pared to the case of say a monolayer of segregation, in 
which the local volumes of the impurities are much more 
relaxed. However, it does not seem possible or useful in 
the step from C to B to try to separate this from the 
chemical contribution, which contains the energy of the 
impurity-host bonds that are broken. Note that our IS 
and CC mechanisms are respectively the "me chanic al" 
and "chemical" contributions identified in refs. m 
the context of interstitial impurities. 

What would be the likely effect of the contributions 
just considered on the work of separation? The HR mech- 
anism should decrease Wscp since removal of a host atom 
cuts more bonds at a grain boundary than at a surface. 
Conversely, the CC step is expected to increase Wsop, 
although the necessity to squeeze an oversized impurity 
while inserting it into a grain boundary acts so as to 
reduce the effect. The reduction should be most pro- 
nounced for tight grain boundary sites or large and soft 
impurity species. (By "large" and "soft" we imply a mea- 
sure such as found in table IVTl below.) Finally, the SS 
step must decrease Wscp for an oversized impurity, as the 
only way for such an impurity to accommodate itself at 
the grain boundary is to push the grains apart, whereas 
at a free surface it may protrude into the vacuum without 
significantly distorting its surrounding host atoms. 

The interstitial impurity case reduces to the compe- 
tition between only two mechanisms (IS and CC); the 
result can be of either signii2i2S42& Interestingly, the IS 
mechanism tends to strengthen the boundary for small 
impurities (H, B, C, N, and O) as they distort free sur- 
faces more than grain boundaries by occupying subsur- 
face sites, whereas the effect of the CC mechanism can 
be either positive or negative depending on the nature 
of chemical bonds formed between impurities and host 
atoms Indeed, in ref "2^ it is found that B and C 
enhance the cohesion of grain boundaries in Mo, while N 
and O reduce the work of separation; the contribution of 
the CC mechanism also acts to lower Wscp if the impurity 
is C, O or N, while only for B does the CC mechanism 
lead to cohesion enhancement at their grain boundary in 
Mo. 

We show in what follows that one can carry the system 
around the gedanken path in figure ^ explicitly by cal- 
culating the quantities defined by equations Q-linjl and 
hence identify the magnitudes of the "size" and "chemi- 
cal" contributions to grain boundary decohesion. 
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III. CALCULATION DETAILS 

We have calculated the atomic and electronic structure 
of the E5 grain boundary in copper from first principles 
using the full potential LMTO method as implemented in 
the NFP codci^ In this code, wavefunctions and charge 
density are expanded in smooth Hankel functions hav- 
ing an energy —k^ and a smoothing radius Tsm defined 
such that beyond smooth and singular Hankel func- 
tions are not appreciably different. For each atom type 
we employ two such envelope functions for each angular 
momentum quantum number up to some £max- These 
are augmented inside atomic spheres with numerical so- 
lutions of the radial Schrodinger equation; all relativistic 
effects except spin orbit coupling are included, the core 
is permitted to relax to self consistency but is treated as 
spherical. In the interstitial region outside the atomic 
spheres, the Poisson equation is solved analytically in 
the smooth Hankel basis and the exchange-correlation 
energy and potential are obtained through a real space 
grid representation. We use the local density approxi- 
mation (LDA) in the parameterisation of von Earth and 
Hedin, modified by Moruzzi et alm^^ Basis function pa- 
rameters and atomic sphere radii used in our study, are 
shown in table ^ For Bi we additionally included a lo- 
calised orbital to allow the Bi 5d semi core to be treated 
as a valence state in addition to the 6d valence state. 
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FIG. 2: Energy-volume curves for bulk Bi. A fragment of the 
respective enthalpy-pressure plot near the bcc-A7 crossover 
is shown in the inset. The energy and enthalpy are measured 
relative to those of Bi(A7) at zero pressure. 

Some bulk properties of Cu and Bi obtained in our 
LDA calculations are listed in table ^1 where they are 
compared with experiment and other LDA calculations. 
Figure 121 shows calculated energy- volume curves for bulk 
Bi. In agreement with experiment, the rhombohedral 



A7 structure has the lowest energy. With pressure in- 
creasing, the AT curve gradually merges into the simple 
cubic curve, the rhombohedral angle steadily increasing 
towards 60° and the internal coordinate u tending to its 
symmetrical value of 0.25. However, at a pressure of 
4.2 GPa the AT curve is intersected by the bcc curve. 
In experiment, the bcc structure becomes stable above 
T.7 GPa, although there are two more phases, Bi-II and 
Bi-III, which appear between AT and bcci^ These are 
rather complex and were not considered in our study. 

The grain boundary was represented in our calcula- 
tions with a periodic supercell containing two identical 
grain boundaries separated by nineteen (310) layers (38 
atoms). The equilibrium width of the grain boundary 
was determined by stretching the supercell along the in- 
terface normal and finding the minimum of the total en- 
ergy of the resulting relaxed supercell as a function of 
the stretch. The lateral dimensions of the supercell were 
fixed by the theoretical equilibrium lattice parameter of 
the bulk fee copper (table Atomic positions were 
fully relaxed in each of the supercells considered until the 
forces dropped below 1 mRy/Bohr. No finite tempera- 
ture effects were taken into account. As the equilibrium 
structure of the S5(310) [001] grain boundary in Cu is well 
known from previous theoretical^ and experimental^^ 
studies, in our calculations we have not attempted to fur- 
ther optimise the atomic structure of the grain boundary 
by shifting the adjacent grains laterally. To minimise 
systematic errors in calculating the work of separation, 
the energy of free surfaces was obtained using the same 
supercell, in which about half of the copper layers were 
replaced with the vacuum so that no grain boundary re- 
mains. For the same reason, we filled the whole super- 
cell with atoms (so called "space filling slab" geometry) 
for the reference bulk calculations. For calculating grain 
boundary and surface excess free energies, equation |(2J), 
at different coverages and hence different stoichiometrics 
of the slabs we employ standard formulaa^i*^ taking the 
chemical potentials of the components from bulk calcu- 
lations. 

Calculations were done for impurity atoms (Bi, Na, or 
Ag) replacing 0, 1 ("loose" site, see figure|2)l, and 2 atoms 
at the grain boundary plane, which we refer to as 0, 0.5 
ML, and 1 ML coverage, respectively. In addition, we 
considered the grain boundary with Bi at 0.25 ML cov- 
erage by doubling the supercell along the [001] tilt axis 
and replacing every other Cu atom in the "loose" posi- 
tion with Bi. In all supercell calculations, a 16x10x2 
Monkhorst-Pack mesh of fc-points was employed with a 
Methfessel-Paxton broadening factor^^ of 5 niRy. To 
represent the smooth interstitial electronic density^^i a 
24x38x152 real space mesh was employed. In practi- 
cal terms, a significant speed up of the calculations was 
achieved by performing the relaxation first with a some- 
what reduced set of parameters (8x4x2 fc-point mesh 
and a 18x25x96 real space mesh) employing the crite- 
rion that forces were stabilised within 1-2 mRy/Bohr, 
and then further optimising the positions using the finer 
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TABLE I: Double basis set parameters used in tlie present study: atomic radius, Va, smoothing radii Vsm., and energetic 
parameters — for each angular momentum £. Radii are chosen so that atomic spheres would not overlap by more than 
few percent in all configurations considered. 



Species 






r 










, Ry 






Bohr 


i^O 


i^l 




^ = 3 


^ = 


I = 1 


i = 2 


^ = 3 


Cu 


2.238 


0.761 


0.667 


0.442 




-0.600 


-0.400 


-0.400 








0.761 


0.831 


0.417 




-1.000 


-1.000 


-1.000 




Bi 


2.450 


0.801 


0.756 


0.486 


0.800 


-0.911 


-0.374 


-0.400 


-0.400 






0.801 


0.756 


0.486 


0.800 


-3.000 




-3.000 


-1.000 


Na 


2.450 


0.569 


0.582 


0.469 




-1.636 


-0.400 


-0.400 








0.569 


0.582 


0.469 




-2.000 


-2.000 


-2.000 




Ag 


2.450 


0.528 


0.551 


0.408 




-1.031 


-0.400 


-0.555 








0.528 


0.551 


0.408 




-3.000 


-3.000 


-3.000 





TABLE II: Bulk properties of Cu and Bi: lattice constant, oo, rhombohedral angle, a, and internal coordinate, u, of the 
equilibrium crystal structure (fee for Cu, and A7 for Bi), structural energy differences, A_E/3 with respect to that structure, 
bulk modulus, B, and elastic constants, C' = (cn — Ci2)/2 and C44, obtained in the present study, other scalar relativistic 
(unless indicated otherwise) LDA calculations, and experiment. 







Cu 






Bi 




our LDA 


other LDA 


Exp. 


our LDA 


other LDA 


Exp. 


ao, A 


3.5270 


3.53", 3.52\ 


3.6146^° 


4.6211 


4.62", 4.528"= 


4.7460^ 






3.58^ 3.52^* 










a 








58.25° 


58.23°", 58.93°*= 


57.23°i^ 


u 








0.2370 


0.238", 0.238*= 


0.234i^ 


AEfcc, mRy 











11.4 






AEbcc, mRy 


2.7 


2.9", 3.2^' 




7.0 






AEsc, mRy 


39.4 


40.0-'^ 




1.5 






B, GPa 


185.9 


190", 192^ 




43.4 


46", 51" 


38.3^ 






153^ 190'* 










C, GPa 


31.2 


35", 27.2^ 25'* 


23.5S 








C44, GPa 


98.9 


99", 86^ 80'' 


7533 









" Ultra soft pseudopotenti alSj Refs. |l3l34| . 
' LAPW, ref. Ha. 
Full potential LMTO, non-relativistic, ref. |36 

** Full potential LAPW, ref. Ha. 
*= Norm conserving pseudopotentials, ref. |3^ . 
^ LAPW, reL i39|l. 



set of meshes. A typical difference between the work of 
separation and grain boundary expansion obtained with 
either set of parameters, 0.05 J/m^ and 2 pm, respec- 
tively, can be taken as a "pessimistic" estimation of the 
convergence of our calculations with respect to the num- 
ber of fc-points and to the size of the real space mesh. 

Fully relaxed calculations of the impurity atoms in Cu 
bulk were done using 32 atom cubic cells with an 8x8x8 
A;-point mesh and a 48x48x48 real space mesh. 



IV. PURE COPPER E5 GRAIN BOUNDARY 

The atomic structure of the relaxed S5(310)[001] sym- 
metric tilt grain boundary in pure Cu obtained in the 



present study, is shown in figure as viewed along the 
[001] direction. It contains two alternating (001) planes 
separated by half the lattice constant, oq, shown with 
larger and smaller circles. It can be imagined as a mirror 
related twin boundary structure in which one of the two 
(310) boundary planes in either grain bordering the mir- 
ror plane is removed, while the other one collapses into 
the grain boundary plane. The same structure can be 
obtained by arranging a suitable lateral shift for the up- 
per half of the original twin with subsequent relaxations^ 
One ends up with two non-equivalent atoms at the grain 
boundary plane, one of which being in the "loose" en- 
vironment and the other being in the "tight" environ- 
ment. The former provides a natural substitutional site 
for an oversized impurity. A recent experimental STEM 
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FIG. 3: Relaxed atomic structure of the Cu E5(310)[001] grain boundary: (a) Pure copper; (b) with 0.5 ML of Bi; and (c) 
with 1 ML of Bi at the grain boundary plane. Notice that in the pure grain boundary (a) there are two inequivalent atomic 
positions in the boundary plane. We refer to these as "loose" (the larger circle) and "tight" (the smaller circle) to illustrate the 
average distance of the nearest neighbours. Naturally we expect a larger atom preferentially to select the loose site to occupy. 



TABLE 111: The energy of the (310) surface, 7^, and of the 
E5(310)[001] grain boundary, 7^6, the grain boundary width, 
w (grain boundary excess volume per unit area) , and the work 
of separation, Wsep, in pure Cu. Experimental data and the 
results of other calculations are given in parentheses. For 
comparison, the results for the E19a(331)[110] grain boundary 
and (331) surface— are also shown. 



Impurity excess (ML) 
0.25 0.5 0.75 



Pure Cu 


E5 


ElOa'' 


w, A 


0.28 


0.22 


7s, J/m^ 


2.21 (1.77") 


2.07 


796, J/m^ 


1.07 (0.980\0.8S 


8"=) 1.04 




3.35 


3.10 



" Exp. polycryst., ref 44]. 
^ Finnis-Sinclair many-body potentials, ref. (45l| . 
EAM, ref. US- Ultra soft pseudopotentials, ref. O. 



image of the grain boundary^^ is very similar to that 
depicted in figure The calculated grain boundary en- 
ergy (table HTHl is in good agreement with semi empirical 
data. In addition, the S5 results listed in the table are 
close to those of ab initio calculations obtained for the 
E19a(331)[110] grain boundary in pure Cu.^"^ The works 
of separation of both grain boundaries are sufficiently 
large to ensure that the boundaries are ductile according 
to the criterion Wsep > ^disi — 1~2 J/m^ as they should 
be. 



V. COPPER E5 GRAIN BOUNDARY WITH 
BISMUTH 

The relaxed atomic positions of the E5 grain bound- 
ary with 0.5 ML of Bi are shown in figure those for 
the 0.25 ML of Bi are very similar except that the grain 
boundary expansion is about half as large (see table lIVp . 

We define the grain boundary expansion as the differ- 
ence per grain boundary between the lengths of supercell 
with and without impurity atoms, noting that supercells 
contain the same number of atoms since our impurities 
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FIG. 4: Dependence of the work of separation on the impu- 
rity excess at the E5(310)[001] Cu grain boundary: Bi (open 
circles), Na (filled circles), and Ag (filled squares). 



are substitutional. In view of the large discrepancy in 
the grain boundary expansion of 50.9 pm for the 0.5 ML 
Bi segregated grain boundary obtained in the present 
study compared with the calculation of Duscher et al.^ 
namely 7-10 pm, using the ultra soft pseudopotentials 
(USPP) within the VASP code,''^ we have additionally 
computed the expansion of one of the E19a boundaries 
with 1 ML of Bi, for which we already knew the USPP 
result obtained with the VASP code.^"* As the NFP and 
VASP expansions were found to agree within 6% (102 pm 
and 108 pm, respectively) we can safely conclude that the 
grain boundary expansions listed in table HVI are the cor- 
rect results within the local density approximation. 

The work of separation and the grain boundary expan- 
sion in table HVl are also plotted in figs. ^ and |S1 against 
the Bi excess at the boundary. Up to 0.5 ML of Bi, both 
quantities comfortably fall on a straight line, meaning 



9 



TABLE IV: Grain boundary expansion and work of separation of the E5(310)[001] Cu grain boundary with impurities. Grain 
boundary coverage of 0.5 ML and surface coverage of 1 ML correspond to an areal density of 5.08 impurity atoms per nm^. 



Cu Bi Ag Na 

Impurity excess, ML - 0.25 0.5 1 0.5 1 0.5 1 

GB expansion 5, pm 25.7 50.9 159.2 13.0 44.5 31.1 154.6 
Segregation energy, eV: 

to the (310) surface 3.11 3.01 0.77 0.77 3.00 2.70 

to the grain boundary 1.69 1.63 1.54 0.79 0.46 1.34 1.26 
W'scp, J/m^: 

into 0.5 + 0.5 surfaces 3.35 2.15 0.96 3.37 2.84 2.00 1.00 

into 1 + surfaces 3.35 2.77 2.23 1.85 3.36 2.24 



Impurity excess (ML) 
0.25 0.5 0.75 1 
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FIG. 5: Same as figure |3 but for the grain boundary expan- 
sion. 



that while Bi atoms fill the "loose" sites their effect is 
additive. 

Although the work of separation decreases as Bi pro- 
gressively substitutes Cu at the "loose" sites of the grain 
boundary, even at 0.5 ML it remains larger than Qdisi, 
which implies that up to that coverage, the boundary is 
ductile. This will have been the situation in the work of 
Duscher et al.^^ From experiment it is known that the 
S5 grain boundary is severely embrittled by Bii-^-'^'' We 
therefore proceeded to add more Bi at the grain bound- 
ary by replacing Cu both at the "loose" and "tight" sites 
(1 ML). In this calculation we did not fully optimise the 
grain boundary structure with respect to rigid lateral 
shifts of the grains, however we removed the mirror sym- 
metry constraint so that the grains could move during the 
atomic relaxation if they wished. As can be seen in fig- 
ure the relaxed structure of the grain boundary does 
show a significant lateral displacement of the grains; in 
addition the Bi layer substantially buckles giving rise to 
large increase of the grain boundary width. The buckling 



of the grain boundary plane is due to large Bi atoms be- 
ing too close to each other, perhaps the structure shown 
in figure |3t can be considered as a precursor of the grain 
boundary faceting which is often detected at Cu grain 
boundaries with Bi^Si^SiS Specifically in the E5 grain 
boundary, some "roughening" has been reported^^ 

The work of separation of the Cu grain boundary with 
1 ML of Bi drops below 1 J/m^ which we take as an 
indication that the threshold concentration of the impu- 
rity required to embrittle the grain boundary, has been 
reached. Note in passing that the way in which impurity 
atoms distribute themselves between the surfaces after 
cleavage turns out to be of little importance; the case of 
equal distribution (denoted as "0.5 -I- 0.5" in table llV|l . 
is marginally lower in energy than the case when all Bi 
atoms reside on one of the surfaces ("1 + 0" in table llVji . 
The exception is the "1 -I- 0" separation of the 1 ML grain 
boundary demonstrating a sizeable increase of Wsop- 

The segregation energies of Bi to the grain boundary 
depend weakly on Bi concentration there (which is just 
another way to say that the Wsop(-r') dependence is lin- 
ear) and are significantly smaller than the segregation 
energies to the free (310) surface, as is expected for an 
embrittling impurity (see Sec. P). Our grain boundary 
segregation energies are comparable to those obtained 
for the SI 9a boundary^ namely 1.32 eV for 0.5 ML 
and 1.58 eV for 1 ML per impurity atom. The experi- 
mentally reported Gibbs segregation free energy to the 
S5 boundary at T = 700°C is Gseg = 0.92±0.11 eV— 
markedly smaller than our theoretical values. This was, 
however, the largest segregation energy among the all 
grain boundaries studied in ref. (45i | . The difference may 
well be due to the effect of vibrational entropy. McLean^ 
argued that the segregation to grain boundaries of so- 
lute atoms sufficiently larger or smaller than the host 
atoms, results in the atomic packing becoming denser at 
the grain boundary, hence the vibrational entropy of the 
grain boundary region may decrease. Decrease of the 
vibrational entropy has been recently predicted for Ag 
segregating to the S5 grain boundary plane in CuS 
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TABLE V: Enthalpy of solution, Hs, and the relative dilation 
volume, nd/f^o, per impurity atom of Bi, Ag, and Na in bulk 
fee Cu. Sid is defined as the difference of the equilibrium 
volumes of the cell with and without impurity at constant 
number of atoms, and f2o is the equilibrium atomic volume 
of pure Cu. The last column gives the formation enthalpy (in 
eV) and formation volume (in units of Q.q) of a vacancy in 
bulk Cu. 



Impurity: 



Bi 



Ag 



Na 



Vac 



Ha, eV: this work 1.66 

other calc. 1.7" 
exp. 

f2ti/r2o: this work 1.40 



0.50 
0.83* 

o.3r 

0.54 0.67 
other calc. 1.3" 0.64^ 0.52'' 
exp. 0.45^ O.SS'' 



1.73 1.27 

1.33''-'', 1.29^ 
1.28'' 
0.64 
0.70*'''' 
0.75'^ 



" Ultra soft DseudoDotentials. R.efs. |l3l34| . 
LSGF, unrelaxecLref. IH. 
Exp., ref. IHIL 
FP-KKR, refTHl. 
Exp., ref. 113. f FP-LMTO, ref. El- 
FP-KKR, ref. 
^ Exp., ref. Eg 



VI. COPPER E5 GRAIN BOUNDARY WITH 
SODIUM AND SILVER 



If Bi embrittlement is a purely structural effect which 
has nothing to do with Bi's chemical identity, then the 
same behaviour is expected from other elements which, 
say, are of about the same size as Bi and are not engaged 
in any strong chemical bonding with the host atoms (see 
section|n]for a more detailed discussion of structural em- 
brittlement). In ref. it was shown that Pb behaves 
very similarly to Bi. One can argue, however, that since 
Bi and Pb are neighbours in the Periodic Table, such 
similarity is indeed expected. We therefore decided to 
pick an element as different to Bi as possible, from which 
only structural embrittlement is expected. For this rea- 
son, we chose the alkali metal Na. Like Bi it does not 
form any stable compounds with Cu and has vanishingly 
small solubility. Indeed, the calculated enthalpy of so- 
lution of Na in fee Cu is as large as that of Bi (table lv)l . 
In addition, the lattice constants of Na and Bi, both cal- 
culated in a hypothetical fee structure, are close to each 
other and much larger than the lattice constant of Cu 
(table lyijl . Na, however, is a much softer species than Bi 
as is exemplified by its smaller bulk modulus (table rvT|) 
and solute dilation volume (table W\ . 

The other element that we have chosen for compari- 
son is Ag for it is known that silver does not embrittle 
Cu. The lattice constant of fee Ag is nevertheless larger 
than that of Cu, and the heat of solution of Ag in Cu is 
positive, although these are not as extreme as for either 
Bi or Na. The experimental Ag-Cu phase diagram shows 
noticeable solubility of Ag in Cu (up to 5%), but without 
the formation of stable compoundsiS 



TABLE VI: Lattice constant, ao, and bulk modulus, B, of the 
impurity species calculated in bulk fee structure. These data 
may be used as a measure of the relative "size" and "softness" 
of atoms. 



Impurity: 


Cu Bi Ag 


Na 


ao, A 


3.5270 4.9048 4.0218 


4.9779 


B, GPa 


185.9 64.7 135.3 


15.8 



TABLE VII: Works of separation H4cp (in J/m^) correspond- 
ing to configurations A-D in figure for impurity coverages 
of 0.5 and 1 ML. Configuration A (pure boundary) of course 
does not depend on impurity coverage. 



Impurity 




0.5 ML 


1 ML 




A 


BCD 


B 


C 


D 


Bi 


3.35 


2.15 1.81 2.20 


0.96 


0.13 


1.38 


Na 


3.35 


2.00 1.97 2.20 


1.00 


0.14 


1.38 


Ag 


3.35 


3.37 2.17 2.20 


2.84 


1.01 


1.38 



The grain boundary calculations with Na and Ag were 
done in the same way as for Bi. The expansion of the 
grain boundary with 0.5 ML of Na is smaller compared 
to Bi, however the grain boundary with 1.0 ML of Na 
expands by almost the same amount as with Bi (see ta- 
ble ^3 a-nd figure IHl). This can be understood assuming 
that Na is a softer species than Bi, but of slightly larger 
size (see table IVl)l and taking into account that the fall 
off of the mutual attraction between Cu grains with dis- 
tance is faster than linear. In all other respects Bi and Na 
behave strikingly similarly (see table lIVl and figure^)). In 
particular, we observe that Na renders the grain bound- 
ary brittle at 1 ML coverage. This is our prediction, as no 
experimental studies of Na at Cu grain boundaries seem 
to exist. However, there are indications that Na can em- 
brittle Cu in the context of liquid metal embrittlement iSi 

The similarity between Bi and Na provides evidence 
that charge transfer effects are of negligible importance. 
If for example the major cause of the embrittlement were 
weakening of Cu-Cu bonds due to transfer of electrons 
from Bi to Cu as proposed in ref. then a free elec- 
tron metal such as Na would have led to even stronger 
embrittlement, say, at 1 ML coverage where the grain 
boundary expansion for Bi and Na is the same. This is 
not observed in our calculations. 

Silver, on the other hand, does not bring about any 
significant decrease of Wscp even at 1 ML (table IIVI and 
figure 01). In fact, at 0.5 ML coverage it even slightly 
strengthens the boundary, despite a small expansion of 
the latter. Such a peculiar situation arises because Ag 
atoms turn out to be of just the right size to replace 
Cu at the "loose" sites. The resulting denser packing at 
the grain boundary thus compensates the replacement of 
Cu-Cu bonds with slightly weaker Cu-Ag bonds. Inter- 
estingly the energy of Ag segregation to the grain bound- 
ary decreases in going from 0.5 ML to 1 ML: this simply 
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TABLE VIII: Impurity segregation energies (in eV per impurity atom) related to W^sep in table IVTTl by virtue of equation 
The sign convention is such that a positive energy means that the impurity wants to segregate. Note that the impurity in 
configurations C and D is actually a vacancy, both in the bulk and at the interface (see section [TTl . 



Impurity 0.5 ML 


1 ML 


BCD 


BCD 


surf gb surf gb surf gb 


surf gb surf gb surf gb 


Bi 3.11 1.63 1.18 -0.71 1.20 -0.22 
Na 3.00 1.34 1.24 -0.46 1.20 -0.22 
Ag 0.77 0.79 1.21 -0.25 1.20 -0.22 


3.01 1.54 1.23 -0.75 1.22 0.02 
2.70 1.26 1.27 -0.70 1.22 0.02 
0.77 0.46 1.25 -0.19 1.22 0.02 



TABLE IX: The changes in Wscp (in J/m^) corresponding 
to Substitutional Structure (SS), Host Removal (HR), and 
Chemical + Compressed Impurity (CC) mechanisms as de- 
fined in figure Q Negative numbers indicate a decrease in 
Wsup ■ For Bi and Na, the HR mechanism provides the largest 
contribution. However, at 1 ML its relative importance di- 
minishes with the SS mechanism contributing nearly half of 
the overall decrease of Wscp- The CC mechanism always acts 
so as to increase Wsep. For Ag, the dominant contributions 
come from the HR and CC mechanisms which are nearly com- 
pensating. (This is an almost trivial result. If we employed 
Cu as the "impurity" then HR and CC would both be large 



and would 


exactly cancel, while SS would be zero.) 


Impurity 


0.5 ML 


1 ML 


Total SS HR CC 


Total SS HR CC 


Bi 


-1.20 -0.39 -1.15 0.34 


-2.39 -1.25 -1.97 0.83 


Na 


-1.35 -0.23 -1.15 0.03 


-2.35 -1.25 -1.97 0.86 


Ag 


0.02 -0.04 -1.15 1.20 


-0.51 -0.37 -1.97 1.83 



reflects the attractive nature of the "loose" site to the 
slightly larger Ag atom. 

VII. HOW DOES BISMUTH EMBRITTLE 
COPPER? 

In accordance with the analysis presented in Sec. ^ 
we have created intermediate configurations C and D de- 
fined in Sec. Ill Al which together with fully relaxed pure 
(A) and segregated systems (B) form a complete set of 
data required. The works of separation for each of these 
configurations corresponding to 0.5 ML and 1 ML cov- 
erages are listed in table IVIll These works of separation 
are further split into surface and grain boundary contri- 
butions in terms of the segregation energies in table lVlill 
Finally, table IIXI shows the changes of Wscp invoked by 
transitions between the configurations A-D that give the 
contributions of the Substitutional Structure (SS), Host 
Removal (HR), and Chemical -I- Compressed Impurity 
(CC) mechanisms, see equations Ql-© in Sec. Ill HI 

For the S5 grain boundary containing 0.5 ML of Bi, 
which is not yet embrittled, we observe that removal of 
host atoms from pure boundary sites that will become 
occupied by Bi (the HR mechanism) strongly dominates 



and almost fully accounts for the total decrease of Wgop 
in table 11X1 The effect of mechanical distortion (the 
SS mechanism) is modest and is similar in magnitude 
to the effect of insertion of Bi back into a vacant site — 
the CC mechanism. This mechanism always acts so as to 
strengthen the boundary. This is exactly what one would 
expect from Ag but perhaps not from Bi if one subscribed 
to the "electronic" point of view. Since the Compressed 
Impurity contribution can only decrease Wsep, yet CC 
is positive, it is clear that the Chemical mechanism is a 
large effect tending to increase the work of separation. 
Even though the grain boundary C is pre-streched, the 
fact that replacing a vacancy even with Bi increases the 
intergranular cohesion is yet another demonstration that 
the electronic grain boundary weakening mechanism of 
embrittlement can be safely ruled out. The introduction 
of Bi divorced from any other effect does not weaken the 
boundary — it strengthens it. 

Increasing Bi coverage to 1.0 ML results in a much 
larger expansion of the grain boundary. The SS contri- 
bution increases substantially giving almost half of the 
overall decrease of Wsep- The HR term still dominates, 
since twice as many Cu atoms are removed from the 
grain boundary plane. Correspondingly, the positive CC 
contribution also increases, but the sum of all three is 
large enough to render the grain boundary brittle. We 
therefore conclude that Bi embrittlement of the Cu grain 
boundary may be attributed to both the Host Removal 
and Substitutional Structure mechanisms with the rela- 
tive weight of the latter increasing as the coverage in- 
creases. 

The effect of Na is strikingly similar to that of Bi (ta- 
ble |^J, especially for 1 ML coverage where not only the 
total effect but even the contributions of each mecha- 
nism are almost identical. Hence everything said about 
Bi can be transferred to Na. At 0.5 ML, the fact that 
sodium atoms are "softer" but "larger" (see table IVI|) 
becomes noticeable with some interesting consequences. 
The grain boundary expansion and hence the SS contri- 
bution is smaller but the total decrease of Wscp is larger 
for Na than for Bi. This is because the CC contribu- 
tion is almost zero. It would not be correct to conclude 
that the bonding between Na and Cu is very weak for 
the effect of CC at 1 ML is as large as for Bi. This is 
a perfect demonstration of the effect of the Compressed 
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Impurity mechanism which combines with the Chemical 
contribution having opposite sign and nearly cancels it 
out. 

Silver does not embrittle Cu grain boundaries, and 
indeed we observe that the HR and CC contributions 
strongly dominate but almost cancel each other (but see 
the caption to table llX|) . At 0.5 ML the "loose" site is 
so attractive for a slightly larger Ag atom that it binds 
there more strongly to other host atoms than even Cu 
itself does, with negligible grain boundary expansion. At 
1 ML silver has to replace Cu at "tight" sites. As a result, 
HR + CC changes sign and the grain boundary expands. 
These combine to decrease W^scp by very small amount 
which is insufhcient to lead to embrittlement. 



VIII. SUMMARY AND CONCLUSIONS 

We have investigated atomic level mechanisms in the 
embrittlement of the E5(310)[001] grain boundary in Cu 
by Bi atoms. To this end, we did first principles calcula- 
tions of the grain boundary with up to 1 ML of impurity 
atoms: Bi, Na, and Ag. We had two principal aims. 

1. Our first was to settle the controversy concerning 
whether embrittlement is a structural or an elec- 
tronic effect. We can reconcile our earlier claim in 
ref. 01 for a structural effect in the I]19a grain 
boundary with apparently contradictory results for 
the E5 grain boundary in ref. |l5j by new calcula- 
tions of the latter which serve to reinforce the con- 
clusions reached in ref. Il3|. We believe the grain 
boundary studied in ref. [l5j was not in fact em- 
brittled. 

In our view the notion that metallic bonds may be 
either weakened or strengthened by charge transfer 
is not tenable. Confusion has arisen because calcu- 
lations supporting these arguments were done using 
clusters of atoms^iS the effects of screening by the 
metallic electron gas were thereby excluded ^ii^S^ 
furthermore total energies were not calculated. We 
do not on the other hand reject chemical effects, 
but as we have shown these strengthen not weaken 
a grain boundary, at least in the case of substi- 
tutional impurities and when defined as here with 
respect to insertion into a vacant site as opposed to 
replacement by a host atom. 

2. Our second aim was to unfold the complex charac- 
ter of embrittlement by dividing the phenomenon 



into independent constituents defined so that for 
any "cause of embrittlement" there exists a sat- 
isfactory thought experiment that isolates it. We 
then proceed to assess in terms of our numerical 
data which of these constituents plays the more 
prominent role. 

Our main conclusions can be summarised as follows. 

1. Segregation of Bi to the S5 Cu grain boundary 
leads to a reduction in its work of separation. We 
predict that at least a monolayer will segregate re- 
sulting in a ductile to brittle transition. 

2. We also predict Na to be nearly as good an embrit- 
tler as Bi. A striking similarity between the em- 
brittling propensity of free electron metal Na and 
semi metallic Bi is strong evidence that the embrit- 
tlement is primarily a structural effect. 

3. Ag is shown not to embrittle the grain boundary, 
in agreement with experiment. In fact at 0.5 ML it 
acts as a modest cohesion enhancer at the S5 grain 
boundary. This is really a size effect: the larger Ag 
atom fits well into the "loose site" at the S5 grain 
boundary. 

4. We suggest an unambiguous, if not unique, thought 
experiment applicable to both substitutional and 
interstitial impurities, which allows one to quantify 
contributions to the changes in Wsop coming from 
the following sources: (i) depletion in the interface 
region (surface or grain boundary) of host atoms 
(for substitutional impurities only); (ii) stretching 
of the host-host bonds at the interface; and {Hi) 
insertion of the impurity into an empty site of pre- 
stretched interface. 

We find the first mechanism to be mostly respon- 
sible for the reduction of the work of separation of 
the boundary with 0.5 ML of Bi or Na. At 1 ML the 
boundary becomes brittle due to the combined ac- 
tion of the first two mechanisms, whereas the third 
mechanism acts in the opposite sense. 
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